3 H production via neutron-neutron-deuteron recombination 
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We study the recombination of two neutrons and deuteron into neutron and 3 H using realistic 
nucleon-nucleon potential models. Exact Alt, Grassberger, and Sandhas equations for the four- 
nucleon transition operators are solved in the momentum-space framework using the complex-energy 
method with special integration weights. We find that at astrophysical or laboratory neutron den- 
sities the production of 3 H via the neutron-neutron-deuteron recombination is much slower as com- 
pared to the radiative neutron-deuteron capture. We also calculate neutron- 3 H elastic and total 
cross sections. 
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I. INTRODUCTION 

The nonrelativistic quantum mechanics solution of the 
four-nucleon scattering problem has, in the past five 
years, reached a level of sophistication and numerical ac- 
curacy that makes it a natural theoretical laboratory to 
study nucleon-nucleon (NN) force models with the same 
confidence as one has used the three-nucleon system in 
the past. This has been demonstrated in a recent bench- 
mark for n- 3 H and p- 3 He elastic scattering observables 
[lj, where three different theoretical frameworks have 
been compared, namely, the hyperspherical harmonics 

(HH) expansion method @, Q , the Faddeev-Yakubovsky 
(FY) equations [H for the wave function components in 
coordinate space [U, and the Alt, Grassberger and 
Sandhas (AGS) equations [3] for transition matrices that 
were solved in the momentum space [HQ. All methods 
include not only the hadronic NN interaction, but also 
the Coulomb repulsion between protons. While the first 
two methods have the advantage of being able to deal 
with charged-particle reactions at very low energies and 
include static three-nucleon forces (3NF), the third one 
is the only method so far to make predictions for multi- 
channel reactions such as d + rf — > d + d, d + d — > ?i + 3 He, 
d + d -» p + 3 H, and p + 3 H — > n + 3 He (and the corre- 
sponding inverse r eactions) [ToL [TTI| . 

In a previous publication [l2| a major step was taken 
in extending the AGS calculations above three- and four- 
cluster breakup thresholds. Owing to the complicated 
analytic structure of the four-body kernel above breakup 
threshold the calculations were performed using the com- 
plex energy method [HI, [l4| whose accuracy and practi- 
cal applicability was greatly improved by a special in- 
tegration method (T^. This allowed us to achieve fully 
converged results for n- 3 H elastic scattering with realis- 
tic NN interactions. We note that the FY calculations 
of n- 3 H elastic scattering have been recently extended as 
well to energies above the four-nucleon breakup threshold 

[HI] , however, using a semi- realistic NN potential limited 
to S- waves. 

In the present work we extend the method of Ref. [lH 
to calculate the neutron-neutron-deuteron (nnd) recom- 
bination into n + 3 H and its time-reverse reaction, i.e., 



the three-cluster breakup n + 3 H — > n + n + d. Al- 
though breakup reactions are usually measured in nu- 
clear physics, the recombination has the advantage that 
its rate is finite at threshold where the breakup cross 
section vanishes due to phase-space factors. Further- 
more, n + n + d — > n + 3 H is the only hadronic re- 
combination reaction in the four-nucleon system that 
at threshold is not suppressed by the Coulomb barrier 
(like n + p + d — > p 
n + n + n + p — > n -f 
environment with neutrons and deuterons and, with re- 
spect to the tritium synthesis, it may be competitive to 
the electromagnetic capture reaction n + d —> 7 + 3 H. 
Thus, one may rise the question at what conditions the 
n + n + d — > n + 3 H recombination would dominate over 
the n + d — > 7 + 3 H radiative capture and to what extent 
it is relevant for astrophysical processes. 

In addition, we also present results for the n+ 3 H elastic 
scattering and study the energy dependence of the total 
n + 3 H cross section. 



f 3 H) or Pauli repulsion (like 
3 H). It can take place in any 



II. 4N SCATTERING EQUATIONS 

We use the time-reversal symmetry to relate the nnd 
recombination amplitude to the three-cluster breakup 
amplitude of the initial n- 3 H state, i.e., 



($l|Ti 3 |$3> = ($3|?3l|$l) 



(1) 



Here is the n- 3 H channel state and |<E>3) is the nnd 
channel state. The advantage is that the three-cluster 
breakup amplitude (^E> 3 IT31 is more directly related 
to the AGS transition operators Up a calculated in our 
previous works 0, [T3 | . Since we use the isospin formal- 
ism where the nucleons are treated as identical fermions, 
there are only two distinct two-cluster partitions, namely, 
P, a = 1 corresponds to the 3+1 partition (12,3)4 
whereas (3 = 2 corresponds to the 2 + 2 partition (12) (34). 
For the initial n- 3 H state we need only Upi, i.e., we solve 
the AGS equations for the four-nucleon transition oper- 
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ators 

U\\ = — (Got Go) 1 -p34 — PzAUiGotGoUn 

+ U 2 G tG U 2U (2a) 

Un = (Got Go)' 1 (I - P 34 ) + (1 - Pb)UiG tGoUu. 

(2b) 

The free resolvent with the complex energy parameter 
Z = E + ie and the free Hamiltonian Hq is 



G = (Z- Ho)- 1 
whereas the pair (12) transition matrix 



t = v + vGot 



(3) 



(4) 



is derived from the respective potential v. The 3+1 and 
2+2 subsystem transition operators are obtained from 
the integral equations 



U a — P a G 1 + P a t Go U a . 



(5) 



the total angular momentum of the (123) subsystem, s y 
and s z are the spins of nucleons 3 and 4, S x , S y , and S z 
are the channel spins of two-, three-, and four-particle 
systems, and J is the total angular momentum with 
the projection M.. We include a large number of 
four-nucleon partial waves, l x ,l y ,l z , jx, jy, Jy < 4 and 
J < 5, such that the results are well converged. The 
complex-energy method [l3| with special integration 
weights [IH is used to treat the singularities of the 
AGS equations ©. To obtain accurate results for the 
breakup amplitude ($3|T3i|$i) near the nnd threshold 
we had to use 0.1 MeV < e < 0.4 MeV that are smaller 
than 1.0 MeV < e < 2.0 MeV used in the elastic 
scattering calculations of Ref. [l2|. However, the need 
for relatively small e values caused no technical problems 
since the integration with special weights [l2| provides 
very accurate treatment of the 3 H pole whereas the 
quasi-singularities due to deuteron pole are located 
in a very narrow region with very small weight, such 
that about 30 grid points for the discretization of each 
momentum variable were sufficient. 



The basis states are antisymmetric under exchange of 
the two nucleons (12). In the 2 + 2 partition the basis 
states have to be antisymmetric also under exchange of 
the two nucleons (34) . The full antisymmetry as required 
for the four-nucleon system is ensured by the permutation 
operators P a b of nucleons a and b with Pi = P\ 2 P 23 + 
P 13 P 23 and P 2 = P 13 P 2i . 

The n- 3 H elastic and inelastic reaction amplitudes at 
the available energy E = e\ + p\/2p,\ are obtained in the 
limit £ — » +0. Here t\ is the 3 H ground state energy, 
Pi is the relative n- 3 H momentum, and p\ = 3mjv/4, 
rriN being the nucleon mass. The elastic scattering am- 
plitude is calculated in Refs. [1, [l2j]. The amplitude for 
the nnd breakup is obtained by the antisymmetrization 
of the general three-cluster breakup amplitude [l(|, re- 
sulting 



(*3|2k|*i) = y3($ 3 |[(l - P 3 4)C/iG tG W 1 
+ U 2 G Q tGoU 21 ]\fa). 



(6) 



Here \<pi) is the Faddeev component of the n- 3 H channel 
state |$i) = (1 + Pi)\(f>i); e% and \<j>\) are obtained by 
solving the bound-state Faddeev equation 



\<h) = G tPM 



(7) 



at £ -t +0. 

We solve the AGS equations ([2]) in the momentum- 
space partial-wave framework. The momentum 
and angular momentum part of the basis states 
are \k x k y k z [l z ({l y [(l x S x )j x s y ]S v }J y s z )S z ] JM) 

for the 3 + 1 configuration and 

\k x ky k z (l z {(l x S X )j X \ly (SyS Z )Sy]jy}S Z )jM) for tllC 

2 + 2. Here k x , k y , and k z are the four-particle Jacobi 
momenta as given in Ref. [l7j . l x , l y , and l z are the 
corresponding orbital angular momenta, j x and j y are 
the total angular momenta of pairs (12) and (34), J y is 



III. RESULTS 

The nnd recombination rate K 3 is defined such that 
the number of recombination events per volume and time 
is K 3 PnPd with p n (pd) being the density of neutrons 
(deuterons). We calculate it as a function of the relative 
nnd kinetic energy E 3 = E — e^, i.e., 



K 3 



(8) 



Here = —2.2246 MeV is the deuteron bound state en- 
ergy, p a y and p a are the reduced masses associated with 
the four-nucleon Jacobi momenta k y and k z . For exam- 
ple, in the 2+2 configuration k y is the relative momentum 
of the two neutrons while k z is the relative momentum 
between the center of mass (cm.) of the two-neutron 
subsystem and the deuteron. The nnd state can be repre- 
sented in both 3+1 and 2+2 configurations equally well; 
liayHa = Tn% /2. The sum in Eq. ([5]) runs over all initial 
and final spin projections m s that are not explicitly indi- 
cated in our notation for the channel states while 93 = 12 
takes care of the spin averaging in the initial nnd state. 
The integral in Eq. ©, up to a factor, determines also 
the total cross section <t 3 for the three-cluster breakup of 
the initial n- 3 H state. Thus, the nnd recombination rate 
can be related to ct 3 as 



g3(^aytla) 3/2 E\ 



■ CT3 



(9) 



where g\ — 4 is the number of 7i- 3 H spin states. Below 
the four-nucleon breakup threshold a 3 can be obtained 
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FIG. 1. (Color online) Differential cross section for elastic 
n- 3 H scattering at 9 MeV neutron energy as a function of 
cm. scattering angle. Results obtained with INOY04 (solid 
curves) and CD Bonn (dashed-dotted curves) potentials are 
compared with the experimental data from Ref. [22| . 



via the optical theorem as a difference between the to- 
tal and elastic cross sections. The equation © cannot 
be used right at the nnd threshold where both E3 and 
(7 3 vanish. For E3 — > the nnd recombination rate (|8]) 
becomes 



4-7T 

K 3 = -(27r)Wi 

93 



where for the channel state 



K*§|r 3 i|*i) 



(10) 



the relative momenta 



0. The most convenient representation for 



is 



a single-component 2+2 state with l y 
i y = and j x = S z = J = 1. 



We study the four-nucleon system using realistic high- 
precision two-nuclcon potentials, namely, the inside- 
nonlocal outside- Yukawa (INOY04) potential by Do- 
leschall @, QJ], the Ar gonne (AV18) potential I19II . the 
charge-dependent Bonn potential (CD Bonn) [20] , and 
its extension CD Bonn + A [2l| allowing for an exci- 
tation of a nucleon to a A isobar and thereby yielding 
effective three- and four-nucleon forces. Among these 
potentials only INOY04 nearly reproduces experimental 
binding energy of 3 H (8.48 MeV), while AV18, CD Bonn 
and CD Bonn + A underbind the 3 H nucleus by 0.86, 
0.48 and 0.20 MeV, respectively. 

First we study the n- 3 H reactions for existing experi- 
mental data. We concentrate on the energy regime rele- 
vant for the nnd recombination, i.e., between the three- 
and four-cluster breakup thresholds. In Fig. [1] we show 
the differential cross section for n- 3 H elastic scattering at 
E n = 9 MeV neutron energy corresponding to E 3 — 0.49 
MeV. The predictions agree well with the experimental 
data of Ref. [22j and are quite insensitive to the choice of 
the potential. Results for n- 3 H elastic scattering above 
the four-cluster breakup threshold up to E n — 22.1 MeV 
are given in Ref. [l2j . 

In Fig. [2] we show the total cross section for n- 3 H scat- 
tering at neutron energies ranging from to 22 MeV 
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FIG. 2. (Color online) Total cross section for n- H scatter- 
ing as a function of the neutron lab energy. Results obtained 
with INOY04 (solid curves) and CD Bonn (dashed-dotted 
curves) potentials are compared with the experimental data 
from Refs. 0,0. 



and compare it to the data of Refs. [23j, |24[ . The three- 
cluster (four-cluster) breakup threshold corresponds to 
E n = 8.35 (11.31) MeV. As already found in Refs. 



[Cam, 



the total n- H cross section around the low- 
energy peak is underpredicted by the traditional two- 
nucleon potentials while the low-momentum or chiral ef- 
fective field theory potentials may provide a better de- 
scription [2^, [27} . Although with increasing energy the 
predictions approach the experimental data, as already 
mentioned [8j, the elastic and total cross section data 
may be inconsistent. In the low-energy regime where the 
inelastic cross section should vanish for E n < 8.35 MeV 
and remain very small at E n = 9 MeV, there is in general 
a better agreement between predictions and experiment 
for the elastic differential cross section than for the to- 
tal cross section which is significantly underestimated by 
theory. A solution to this discrepancy may require new 
measurements in this energy regime. 

In Fig. [3] we study the energy-dependence of the nnd 
recombination rate in the standard form N\K$ where 
Na is the Avogadro's number. We show only INOY04 
predictions as it is the only used potential with correct 
ei and pi values. The results at E 3 — are obtained 
from Eq. (JTUJ) while at E3 > it was more convenient to 
use Eq. ([9|) where 0-3 was calculated using optical theo- 
rem. Thus, for E3 > \ed\ our predictions in Fig. [3] esti- 
mate the upper limit for N^K^ since they assume that 
the four-cluster breakup cross section is much smaller 
than the three-cluster breakup cross section. In the rele- 
vant energy regime < £3 < |ed| the nnd recombination 
rate increases with increasing energy E3 nearly linearly 
due to the increasing contributions of partial waves with 
nonzero orbital angular momentum l z . The threshold 
values N\K® referring to all employed potentials are col- 
lected in Table HI they increase with 3 H binding energy. 
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FIG. 3. (Color online) nnd recombination rate Kz as a 
function of relative kinetic nnd energy E$. Predictions are 
obtained using the INOY04 potential. 
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TABLE I. nnd recombination rate at threshold calculated 
with different two-nucleon potentials. The values for 3 H bind- 
ing energy are listed as well. 



Finally we compare the relative importance of the 
nnd recombination and nd radiative capture. For the 
latter the number of events, i.e., the number of produced 
3 H nuclei per volume and time is K^PnPd with Ki 
being the nd capture rate. The threshold value for it 
given in Ref. Q is N A K% = 66.2cm 3 s" 1 mor 1 . The 
critical density of neutrons at which both processes 
yield comparable contributions to the 3 H production 
in the low-energy (low-temperature) limit is given by 
Pn = K ll K l ~ 2.6 x 10 30 cur 3 . This corresponds to 
the mass density of 4.4 x 10 6 g/cm 3 . Thus, one may 
conclude that at the neutron density available in the 
laboratories (such as National Ignition Facility with 
expected p n ~ 10 22 to 10 25 cm -3 [29[) the nnd recombi- 



nation is entirely irrelevant as well as for the big-bang 
nucleosynthesis where the estimated baryon density is 
even lower. On the other hand, the neutron density in 
core-collapse supernova or neutron stars may be higher 
than p c n by several orders of magnitude but the absence 
of deuterons renders n + d and n + n + d reactions 
irrelevant. However, based on our results one may 
conjecture that at such high densities the three-cluster 
recombination of two neutrons and a heavier nucleus A, 
i.e., n + n + A — > n + {An) might be as important as 
the corresponding radiative capture n + A — > 7 + (An). 
For example, the above reactions with A being 20 Ne are 
relevant for the neon-burning process. 

IV. SUMMARY 

We have solved the four-nucleon AGS equations in the 
energy regime above the three-cluster threshold and stud- 
ied the rate for the recombination reaction n + n + d — > 
n + 3 H. The obtained results show that the nnd recom- 
bination is not competitive with the radiative nd capture 
for the production of tritium at neutron densities avail- 
able in laboratory induced fusion or astrophysical pro- 
cesses. Thus, one may conjecture with a confidence that 
other nucleon-nucleon-deuteron recombination reactions 
(for example, p+p + d — >• p + 3 He that could contribute to 
the hydrogen burning process in stars) , being in addition 
suppressed by the Coulomb repulsion, are inferior to the 
respective nucleon-deuteron radiative capture reactions 
at realistic densities, and that four-nucleon recombina- 
tion reactions are even far less relevant. 

In addition, wc presented results for the n- 3 H elastic 
differential cross section at E n = 9 MeV and the n- 3 H 
total cross section up to E n = 22 MeV. While the elastic 
cross section agrees fairly well with the data, there is a 
disagreement for the total cross section, especially in the 
low-energy regime. This indicates a possible inconsisten- 
cies between the n- 3 H elastic and total cross section data 
and calls for new measurements. 
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